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HEART BEAT SIGNAL ANALYSIS 



The present invention relates to methods and 



5 apparatus for analysing quasi -repetitive signals 

■ 

comprising a series of similar complexes, and in 
particular, but not exclusively, to such methods and 
apparatus suitable for analysing heart beat signals such 
as electrocardiograph signals. 



An electrocardiogram (ECG) is a graph of the voltage 
variations plotted against time. The variations results 
from the depolarisation and depolarisation of the cardiac 
muscle, which produces electrical fields that reach the 

15 surface of the body where electrodes are located Placement 
of multiple electrodes at selected points on the body 
enables a variety of different signals to be picked up 
simultaneously (multiple lead electrocardiography) , from 
which a reasonably complete picture of the electrical 

20 activity of the heart, and any related clinical 

abnormalities can be deduced. A heart beat complex from a 
typical ECG signal is shown in simplified form in figure 
1. The signal comprises many such complexes in series, 
each corresponding to one heart beat. An ECG signal may 

25 be taken from each lead. Typically, twelve leads are 
used. 

The principal features of an ECG complex have widely 
accepted designations. An initial P-wave is followed by 

30 a larger QRS-wave. A T-wave follows the QRS-wave. Various 
minor variations in the shape of a complex, as well as 
missed and extra beats are expected in a healthy subject. 

Heart abnormalities, diseases, and drugs also affect the 
shape and timing of the features of a complex. Trials of 

35 a drug generally include experiments to determine any such 
changes, from which the effect of a drug on the 
electrophysiological activity of the heart can be deduced. 
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Pre-clinical drug trials generally include 
administering a drug to animals using various dosage 
schemes. During an experiment, which typically lasts 
5 several hours, one or several doses of a drug may be 
administered and ECG recordings are made. It is time 
consuming to manually analyse the large amount of ECG data 
resulting from such experiments, and from other situations 
in which long series of heart beats are recorded. 

10 

Automatic analysis of ECG complexes to identify 
various parameters is well known. Analysis of the mutual 
timings of various features of signals from one or more 
ECG leads, for example, is routinely built into ECG 

15 machines for clinical use. Such shape analysis is usually 
carried out on single selected heart beat complexes. 
Automatic shape analysis of the average of a series of 
complexes is also known, for example from US 5,139,027, 
which describes an apparatus which also automatically 

20 looks for features in an ECG signal such as narrow spikes, 
electrical noise and other artifacts and discards them 
before forming an average of a plurality of complexes. 
The averaging of a plurality of complexes before 
performing shape analysis is advantageous in that it 

25 reduces the computational cost of the expensive shape 

analysis process and improves the signal to noise ratio of 
the analysed signal. 

However, the known schemes for automatic analysis of 
30 ECG signals are not ideal in many situations, such as pre- 
clinical drug trials, because they do not adequately allow 
for natural , variations in shape and timing between the 
members of a group of heart beat complexes to be averaged. 

35 The invention seeks to address the above and other 

problems of the related prior art. 
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Accordingly, the invention provides a method of 
analysing a heart beat signal which comprises a series of 
heart beat complexes, the method comprising the steps of: 

selecting a set of said complexes; 

identifying from the set a subset of complexes each 
of which satisfy a similarity criterion with respect to 
members of the set; 

calculating a representative complex from the members 
of the subset; and 

establishing one or more heart beat signal parameters 
from said representative complex. 

The heart beat signal is preferably an ECG signal, 
but the method may also be applied to other types of heart 
beat signal, such as a signal from a pulse oximeter, a 
thoracic motion sensor or an ultrasound device. Each 
heart beat complex preferably corresponds to a single 
heart beat, or portion thereof. Each selected set of 
complexes preferably includes only genuine complexes, 
excluding artifacts of various sorts familiar to the 
skilled person, such as spikes and electrical noise. The 
complexes in any one set are preferably consecutive within 
said signal. Typically, the set will comprise between 
about two and twenty complexes , and preferably between 
four and ten complexes . 

The heart beat signal parameters may be calculated 
from the representative complex using a variety of 
techniques known in the art. For an ECG signal, 
parameters, such as the PQ, QRS and QT intervals, and 
parameters derived from intervals between features on 
different leads may be calculated. 

Preferably, the step of identifying comprises 
identifying from the set a subset of complexes each of 
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which satisfy a similarity criterion with respect to an 
average of the members of the set, where the average could 
be a mean, an unnormalised sum, or a variety of other 
types of average. Preferably, said step of calculating 
5 comprises calculating the average complex of the members 
of the subset, where again the average could be one of 
several types of average. 

Preferably, the method includes a step of analysing 
said heart beat signal to identify and select genuine, or 
valid heart beat complexes to form the above mentioned 
series, before any of steps mentioned above are carried 
out. The method may also include steps of aligning the 
complexes of the series, the set and/or the subset prior 
to operations such as correlation and averaging. 

Preferably, the similarity criterion comprises a 
threshold calculated using the complexes of said set and 
an average of these complexes. For example, the criterion 
could be based on coefficients of correlation between each 
member of the set and the average of the set, the 
criterion comprising a threshold correlation coefficient. 

Such a threshold coefficient may itself be based on the 
distribution of the calculated correlation coefficients 
for the set. However, the similarity criteria preferably 
also comprises a predefined threshold value, for example a 
fixed correlation coefficient of 0.98, above which the 
criterion is always satisfied. This ensures that no 
complexes fail the similarity test if they are all very 
similar . 

Said step of identifying may further comprise 
aligning the complexes of said set, forming a set- 
average complex from said aligned complexes, comparing 
35 each complex of said set to said set-average complex, 

determining if each result of comparing satisfies 
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said similarity criterion, and forming said subset from 
only those complexes of said set for which the similarity 
criterion is satisfied. Preferably, the step of comparing 
comprises correlating each complex of said set with said 
5 set-average complex. 

If said signal is an electrocardiograph signal, said 
step of aligning the complexes of said set preferably 
comprises the step of aligning said complexes on an R-top 
10 feature of each complex, although other alignment schemes 
could also be used. 

The method steps described above may conveniently be 
implemented using computer software for execution on a 

15 suitable computer system. Such software may conveniently 
be written on one or more computer readable media such as 
CDROMs . The invention also provides a computer system 
arranged to carry out the methods . Such a computer system 
may provide, for example, an extract set element adapted 

20 to select a set of complexes from the series of complexes 
to be analysed, a select subset element adapted to 
identify from the set a subset of complexes each of which 
satisfy a similarity criterion with respect to members of 
the set, and a combiner element adapted to calculate a 

25 representative complex from the members of the subset. 

The invention also provides apparatus for analysing a 
heart beat signal comprising a plurality of heart beat 
complexes, the apparatus comprising: an analysis engine 
adapted to automatically calculate a representative heart 

30 beat complex for each of a plurality of preselected 
intervals of said signal; and 

an editor adapted to enable a user to edit calculation 
and/ or heart beat parameters relating to each 
automatically calculated representative heart beat 
35 complex. Typically, each interval contains a set of heart 
beat complexes . The analysis engine is preferably adapted 
to automatically decide which complexes in the set to 
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include in each representative heart beat complex, and the 
editor is preferably adapted to enable the user to 
subsequently change which complexes are included. 

5 The editor may enable the user to control calculation 

parameters which determine whether a particular complex 
should be included, not included, or included subject to 
an automatic or predefined similarity criterion, in the 
corresponding representative heart beat complex. This 

10 representative complex may be an average of the included 
complexes . 

The editor may further provide the user with a 
graphical display of at least some of the heart beat 

15 complexes within a current one of the intervals of said 
signal, and with a graphical display of the current 
representative heart beat complex for the current 
interval. In this way, the user may be provided with 
facilities to validate the representative heart beat 

20 complex and correct it if necessary. 

For example, the editor may be further adapted to 
display, on the display of the current representative 
heart beat complex, one or more feature markers indicating 

25 the locations of one or more automatically determined 

features of said representative heart beat complex, such 
as P, Q, R, S and T features of an ECG complex. The 
editor is preferably further adapted to enable the user to 
move said markers, for example by means of a drag and drop 

30 action using a pointer device. 

Preferably, the analysis engine is adapted to 
calculate one or more heart beat parameters, especially 
shape parameters, of each representative heart beat 
35 complex, and said editor is adapted to display said heart 
beat parameters and to update said display according to 
changes made by said user. Such heart beat parameters 
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could include, for ECG analysis, the duration of the QT, 
PQ and QRS as well as of other intervals between features 
for a single lead or between leads. 

The apparatus is preferably provided by means of a 
computer program and appropriate data elements, which may 
be provided on one or more removable computer readable 
media, executing on a computer. One preferred embodiment 
is implemented using Microsoft Excel (RTM) spreadsheet 
templates, macros and other appropriate data structures 
and files. 

The method and apparatus may be used in a variety of 
circumstances and applications. In the embodiments set 
out below, pre-clinical drug trials are particularly 
mentioned, but other applications include any analysis of 
heart beat data including clinical drug trials, medical 
check ups, patient monitoring in intensive care, during 
and after surgical procedures, in athletic activity 
monitors and so on. 

The methods and apparatus may also be applied to 
repetitive or quasi-repetitive signals other than heart 
beat signals including to signals of medical origin and 
use and to non-medical signals. 

Embodiments of the invention will now be described, 
by way of example only, with reference to the accompanying 
drawings, of which: 

Figure 1 shows features of a typical ECG signal; 

Figure 2 sets out a scheme for analysing a heart beat 
signal ; 

Figure 3 illustrates how steps 14 and 16 of figure 2 
may be put into effect, in more detail; 

Figure 4 illustrates apparatus and/or computer 
program elements for putting the method steps of figures 2 
and 3 into effect; 

Figure 5 sets out a scheme for analysis and 
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validation of a heart beat signal; and 

Figure 6 shows a graphical user interface for 
assisting in the validation step of figure 5. 

Figure 2 illustrates, schematically, a method of 
analysing a heart beat signal 10, such as an ECG signal, 
embodying the invention. The method comprises steps of 
selecting a set of heart beat complexes from within the 
signal (step 12), identifying a subset of those complexes 
which all- match a quality, similarity or mutual similarity 
criterion (step 14) and calculating an average or 
representative complex of the identified subset (step 16) . 

The average complex may be stored and used to calculate 
one or more heart beat parameters (step 18) . If the 
signal is an ECG signal, then these heart beat parameters 
may be, for example, a QRS interval, a PQ interval, a QT 
interval and so on. 

Steps 14 and 16 of figure 2 may be implemented as 
illustrated in figure 3 . The members of the selected set 
of heart beat complexes 20 are first aligned in step 22. 
The alignment may conveniently be carried out using a 
single dominant feature such as the R-peak of an ECG 
complex, or may involve optimization of alignment of more 
then one feature or region of each complex. After 
alignment, the complexes are then summed or averaged in 
step 24. A correlation coefficient between each aligned 
complex and the averaged complex is then calculated in 
step 26. The set of correlation coefficients so obtained 
is used to calculate a threshold parameter in step 28 . 

The threshold parameter calculated in step 28 is used 
in step 30 to discard one or more of the original 
complexes 20 according to a similarity criterion using the 
parameter. The remaining subset of complexes may again be 
aligned, for example with reference to an R-top, in step 
32, and the subset is then used to calculate an average 
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complex, representative of the original complexes 20, in 
step 34. 

A suitable threshold parameter calculated in step 28 
is the value above which 95% of the correlation 
coefficients would be expected to lie, based on the 
average and standard deviation of their actual values. 
The similarity criterion used in step 30 could then be 
such as to discard complexes having a correlation 
coefficient calculated in step 26 falling below this 
threshold parameter, subject to a maximum threshold 
parameter value of 0.98. 

The method illustrated in figures 2 and 3 may be put 
into effect on any suitable computer equipment, in 
particular using suitable software elements which may be 
provided on one or more computer readable media. Suitable 
elements of such software, or computer apparatus suitably 
programmed, are illustrated in figure 4. The arrangement 
may be controlled and the results of the implemented 
process displayed using the apparatus and user interface 
further described below in relation to figure 6. 

Figure 4 illustrates a memory device 40, such as a 
disk drive or RAM, storing an ECG or other heart beat 
signal. The signal may be preprocessed in various ways by 
preprocess element 42. The preprocess element 42 may 
carry out cleaning and filtering processes on the signal, 
and in particular may identify and label valid heart beat 
complexes. Noise artifacts and badly distorted complexes 
are rejected, for example by comparison against templates 
representing acceptable complex shapes, by checking the 
timing of each complex within the signal, or by other 
techniques known in the prior art . 

An extract element 44 selects one or more sets of 
complexes for further analysis, each set preferably 
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containing a fixed number of consecutive complexes. For 
the purposes of pre-clinical drug trials each set could be 
a series of ten complexes selected at fifteen minute 
intervals throughout an experiment, with extra sets being 
5 selected in each period immediately following a drug 
administration. 

An align element 46 aligns the complexes of each set. 
A calculate criterion element 48 calculates the parameter 

10 or parameters, such as a threshold, for deciding which 

complexes of the set to use further, and which to reject, 
for example as described above in connection with steps 
24, 26 and 28 of figure 3. An identify subset element 50 
uses the parameter or parameters from the calculate 

15 criterion element 48 to select a subset of the complexes, 
as described above in connection with step 30 of figure 3. 

Finally, the selected subset is used by averaging element 
52 to calculate a representative complex, and the analyse 
complex element 54 automatically analyses the average 

20 complex to obtain desired heart beat parameters. At each 
stage of the process, and especially at the final analysis 
stage, results may be written back to the memory 40 in a 
manner so to be suitably associated with the original 
signal data. 

25 

Generally, and particularly for pre-clinical drug 
trials, the heart beat signal analysis scheme described 
above is incorporated into a process such as that set out 
in figure 5. An ECG signal is recorded at step €0. 

30 Corresponding experiment information is input at step 62, 
and used to perform an automatic analysis of the ECG 
signal at step 64. The results of the automatic analysis 
are checked, and if necessary adjusted at validation step 
66, and finally an experiment report is produced at step 

35 68. The experiment information may include conditions of 
the experiment, such as time, date, duration, location, 
subject animal, drug, dosage and so on, as well as 
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parameters expressly constraining the subsequent 
calculation step 64 such as the time points at which 
detailed analysis is required, how that analysis is to be 
carried out, and what result heart beat parameters are to 
be obtained. 

■ 

In one embodiment of the invention, steps 62 to 68 of 
figure 5 are implemented using Microsoft Excel (RTM) 
spreadsheet templates and macros, with ECG signal data 
being imported in a compatible file format. Figure 6 
illustrates a graphical user interface provided as part of 
this embodiment to implement the validation and adjustment 
step 66 of figure 5 following the calculation step 64. 
This interface is displayed by a Microsoft Excel (RTM) 
tool executing on a conventional personal computer, 
providing conventional pointing device, keyboard and 
visual display peripherals. 

The duration of the experiment being validated in the 
interface of figure 6 is shown in the time point area 70. 

The time points (in minutes) at which detailed analysis 
of the ECG signal is required are shown. The first time 
point is minus fifteen minutes, because the zero minute 
point is the time of a first drug administration. A 
current time point is indicated by a box 72 surrounding 
the m -15" value. Completed validation at a time point is 
indicated by colour coding in the time point area 70. 

The ECG parameters being validated are based upon up 
to twelve time-parallel signals from up to twelve ECG 
leads, listed in the left hand column of a validation 
table 76. A current ECG lead is indicated by a box 78 
surrounding the lead II row. Colour coding in the 
validation table 7 6 illustrates which beats of which lead 
have been validated, and which leads do not require 
validation. A protocol editor may be available. In said 
protocol editor the user can predefine time points, leads, 
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parameters and so on that are to be used or calculated in 
the ECG editor. The horizontal axis of the lead table is 
representative of the ten beats (or other number as 
appropriate) to be averaged at the current time point. 

The ECG signal of the beats to be used at the current 
time point for the current lead are plotted in a beat plot 
area 80, as multiple curves overlaid according to an 
automatic alignment process. The ten beats shown in the 
beat plot area 80 of figure 6 have been aligned on their 
R-top features. Which beats are plotted can be controlled 
using a beat plot selector 82, which currently indicates 
that all beats of the set except beats which have been 
rejected on the basis of a similarity test (ie used beats) 
are shown. Other options provided by selector 82 are all 
beats of the set, and single selected beats. 

Below the beat plot area is an average beat plot area 
86 which displays, using a horizontal axis aligned with 
that of the beat plot area 80, an average beat calculated 
according to the method illustrated in figures 2 to 4 . 
Also shown are markers 88 (P,Q,S and T) showing the 
positions of corresponding features of the average beat 
which have been automatically calculated, for example 
using beat shape analysis techniques known in the prior 
art. The markers 88, which could mark various different 
features, can be moved by the user. Heart beat parameters 
derived from the current marker positions are given in a 
heart beat parameter area 90. 

Although the average beat and consequent heart beat 
parameters are automatically calculated for all time 
points according to a schedule included with the 
information provided in step 62 of figure 5, the user of 
the interface of figure 6 can subsequently validate and 
intervene in the beat averaging process by means of the 
controls provided in the average beat control area 94. 
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Each row of this area corresponds to one of the beats at 
the current time point. The "Correl" column shows the 
correlation coefficient of each beat, once aligned, with 
the average of all beats of the set, and the "Correl U» 
> column shows the correlation of each beat, once alxgned, 
with the average of all beats currently used to form the 
average beat shown in the average beat plot area 86. 

The "User" column of the average beat control area 94 
0 enables the user to allow the tool to decide whether to 
include a beat in the average according to the 
predetermined scheme, or to override the scheme and choose 
to include or exclude the beat from the average. This 
facility allows a user to reject a distorted or otherwxse 
13 incorrect or undesirable beat, and to include an otherwxse 
rejected beat if necessary. 

Navigation area 98 provides facilities for moving the 
current validation point between leads, beats, time poxnts 
20 and experiments, for accepting and validating data and 
derived parameters and so on. 

Alternatives and variations to the embodiments 
described above will be apparent to the person skilled xn 
25 the art, while still falling within the scope of the 
appended claims and their equivalents. For example, a 
variety of different programming languages could be used 
to implement the invention, on a variety of different 
hardware platforms including dedicated ECG machines, 
30 portable heart beat monitors, and conventional computer 
work stations. A validation interface could include 
various combinations of graphical elements selected from 
those illustrated in figure 6, as well as other elements. 
The validation interface may also be provided as a 
35 separate entity to any automatic processing means, the 
results from which it is to be used to validate. 
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Various criteria may be used to select which heart 
beat complexes to include when calculating a 
representative complex. In some applications only 
complexes deviating to a large extent from the average of 
5 a set will be rejected, while in other applications only 
near identical complexes will be included. 



